NASA Technical Paper 1799 


LOAN con: RF 
AfWL TECHNlCr* 
KIRTLAND AfB, 


Analysis and Simulation of a 
Magnetic Bearing Suspension 
System for a Laboratory Model 
Annular Momentum Control Device 


Nelson J. Groom, Charles T. Woolley, 
and Suresh M. Joshi 


MARCH 1981 


NASA 
TP 
1799 
c. 1 



Ln - 

^ -- 





TECH LIBRARY KAPB, NM 



TECH UBFtARY KAFB, NM 

NASA Technical Paper 1799 



Analysis and Simulation of a 
Magnetic Bearing Suspension 
System for a Laboratory Model 
Annular Momentum Control Device 


Nelson J. Groom and Charles T. Woolley 
Langley Research Center 
Hampton, Virginia 

Suresh M. Joshi 

Vigyan Research Associates, Inc. 

Hatnpton, Virginia 


IM/\SA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Branch 


IL. 


1981 



SUMMARY 


A linear analysis and the results of a nonlinear simulation of a magnetic 
bearing suspension system which uses permanent magnet flux biasing are presented. 
The magnetic bearing suspension is part of a 4068 N-m-s (3000 Ib-ft-sec) labora- 
tory model annular momentum control device (AMCD) . Included in the simulation 
are rigid body rim dynamics, linear and nonlinear axial actuators, linear radial 
actuators, axial and radial rim warp, and power supply and power driver current 
limits. 


INTRODUCTION 

The basic concept of the annular momentum control device (AMCD) is that 
of a rotating annular rim suspended by noncontacting magnetic bearings and 
powered by a noncontacting linear electromagnetic motor. A detailed discussion 
of the rationale for the AMCD configuration and its potential applications are 
presented in reference 1 . Earth-based energy storage applications of the con- 
cept are discussed in references 2 to 4. 

In order to investigate any potential problems in implementing the AMCD 
concept, a laboratory model AMCD was designed and fabricated under contract. 

This hardware has been delivered and preliminary tests have been performed. 

The appendix presents a brief description of the laboratory model, including 
a table of pertinent physical parameters; reference 5 presents a detailed 
description. Reference 6 presents results of static and low-speed dynamic 
tests, which include spin motor torque characteristics and spin motor and mag- 
netic bearing drag losses. Reference 6 also briefly discusses permanent mag- 
netic flux biasing and rim warp. This paper presents an analysis of the 
laboratory model AMCD magnetic bearing suspension system which addresses the 
subject of permanent magnet flux biasing and rim warp in more detail. A digital 
computer simulation, which is used in the analysis of the laboratory model, is 
described in reference 7. Reference 8 presents the development of an analytical 
model of the laboratory model magnetic bearing actuator. 


SYMBOLS 

Values are given in both SI and U.S. Customary Units. The measurements 
and calculations were made in U.S. Customary Units. 

[a] matrix defined by equation (1 0) 

total axial force acting on AMCD rim produced by axial bearings a, 
b, and c 

Fj^ force produced by a given radial bearing 



^Ra'^b'^Rc forces produced by radial bearings a, b, and c, respectively 
force command for a given radial bearing 
Fri ,Pj^ 2 forces on AMCD rim along axes 7 and 2, respectively 

tFx> = f^Xa Fxb Fxcj^ 

^Xa'^Xb'^Xc forces produced by axial bearings b, and c, respectively 
f^A2 torques on AMCD rim about axes 1 and 2r respectively 
G 3 j^ radial bearing transfer function 

Gp forward loop transfer function of axial bearing loop 

Gpi^ forward loop transfer function of radial bearing loop 

gp gap for a given radial bearing 

^Ra'^Rb'^Rc gaps for radial magnetic bearings a, b, and c, respectively 

^RaC'^RbC^^RcC commands for radial magnetic bearings a, b, and c, 

respectively 

9rC command for a given radial bearing 

{gx^ = £Xa 9xb gxcj"^ 

^Xa'^Xb'^Xc 9^ps for axial magnetic bearings a, b,^ and Cf respectively 
^A03 angular momentum stored in rim about spin axis 

[ll identity matrix 

moment of inertia of rim about transverse axes 
moment of inertia of rim about spin axis 
position gain of a given axial bearing control loop 
Kg equivalent electromagnet gain of a given magnetic bearing 

equivalent permanent magnet stiffness of a given magnetic bearing 
Kp rate gain of a given axial bearing control loop 

[m] matrix defined by equation (23) 

m;^ mass of rim 

^CA7 ' ^CA2^ ^CA3 translations along axes 1, 2, and 3, respectively 
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r 


*^CA1 C' ^CA2C' ^CA3C translation commands along axes 1, 2, and 3, 

respectively 

r^i mean radius of rim 

s Laplace variable 

[t] matrix defined by equation (11) 

weight of rim 

®Al '®A2 rotations about axes 1 and 2, respectively 

0A1C'®A2C rotation commands about axes 1 and 2, respectively 

p damping ratio of axial bearing loop with ideal magnetic bearing 

Pg total system damping ratio 

p damping ratio of axial bearing loop with permanent magnet flux biasing 

(j)g spin rate of rim 

damped natural frequency (eq. 38) 

damped natural frequency of system high-frequency poles 

damped natural frequency of system low-frequency poles 

natural frequency of axial bearing loop with ideal magnetic bearing 

natural frequency of axial bearing loop with permanent magnet flux 
biasing 

Matrix notation: 

[ ] rectangular matrix 

[ ] “^ inverse of rectangular matrix 

[ ] transpose of rectangular matrix 

{ } column vector 

I J row vector 

Dots over symbols denote derivatives with respect to time. 
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EQUATIONS OF MOTION 


The equations of motion used in this paper are for a rigid AMCD rim sus- 
pended by magnetic bearings attached to a rigid fixed base. They are described 
in detail in reference 9. The location of the magnetic bearing suspension 
stations with respect to the AMCD axis system is shown in figure 1 . 

The axial gaps, in terms of rim rotations and translations^ are given by 


9xa 

\ 9xb > = 
9Xc 


(1/372) 

-( 1/3/2) rm 


-(V2)r„, 

-(V2)r^ 





1 


®A1 

1 

< 

®A2 

1 


*^CA3 





0) 


where gjjar Uxb^ 9Xc axial gaps for bearing stations a, b, and c, 
respectively, rj^ is the mean radius of the AMCD rim, 0;^1 and 0;^2 rota- 
tions of the rim about axes 1 and 2, respectively, and r^/va is the translation 
of the rim along axis 3. The radial gaps in terms of rim translations are 


9Ra 


1/2 

1/3/2 

9Rb 

) = 

1/2 

-1/3/2 

9rc 


-1 

0 



( 2 ) 


where URa^ 9Rb^ 9 rc radial gaps for bearing stations a, b, and 

c, respectively, and r^Al and r^^2 translations of the AMCD rim along 

axes 1 and 2, respectively. The torques on the rim due to the axial bearing 
forces are 



(\/V2)rn, 

-0/2)rro 


(\/372) Tjn 

-n/2)r„ 


0 

^m 


^Xa 

(Fxb ) 

^Xc 

V J 


(3) 


where and G;^2 torques on the rim about axes 1 and 2, respectively, 

and F^af ^Xb' ^Xc axial forces produced by bearing stations a, 

b, and c, respectively. The radial forces resolved along axes 1 and 2 are 
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1/2 
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where Fri and Fr 2 are forces on the rim along axes 1 and 2 and 

and are the radial forces produced by bearing stations a, b, and c, 

respectively. The rim rotational and axial translation dynamics are 


• • t 

I A® AT = - Hao3®A2 

(5a) 

^A®A2 = Ga2 + H;^03®A1 

(5b) 

™A’rCA3 = ^A ■ Wa 

(5c) 


where is the rim transverse inertia, m^ is the rim mass, H;^q 3 is the 

rim angular momentum about axis 3 (the spin axis) , and is the weight of 

the rim. 


= ^Xa + ^Xb + Fxc 


( 6 ) 


Finally, the radial translation dynamics become 



MAGNETIC BEARING SUSPENSION SYSTEM 

The laboratory model AMCD magnetic bearing suspension system provides 
active positioning control of the rim in both the axial and radial directions. 
The axial and radial suspension systems are independent and are designed 
separately. 


Axial System 

Assuming three magnetic bearing suspension stations equally spaced around 
the rim and ideal rim inertia distribution, it will be shown subsequently that 
axial motions of the rim in each of the bearing stations are uncoupled at zero 
rim speed. That is, axial motion of the rim in one bearing produces no motion 
in the other two bearings. Consequently, at zero momentum, the axial magnetic 
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bearing control system decouples into three identical independent systems. Thus, 
a single design, using a simplified suspended mass model, can be performed for 
each system. Applying this decoupled design approach to the laboratory model 
AMCD, the closed-loop magnetic bearing control system parameters required to 
produce desired system performance at a given rim angular momentum are obtained 
by analyses similar to those in reference 1 but with respect to a fixed base. 

To illustrate this approach, the axial system characteristic equation will 
be developed with the assumption that the magnetic bearing actuators are perfect 
(i.e., force output is equal to force command) and that the force commands are 
functions of rim position and rate. Under these assumptions, the bearing force 
F as a function of bearing gap g at a given bearing station can be written as 


P = Gpg = (Ka + KRS)g 


( 8 ) 


where is a position gain, Kr is rate gain, Gp is the forward loop 

transfer function, and s the Laplace operator. The rim rotational and axial 
translation equations of motion are given in equations (5) . Taking the Laplace 
transform of these equations and putting them in matrix form results in 


Ias2 

Ha03S 

0 


®A1 


^A1 

-Ha03S 

I AS 2 

0 

< 

®A2 

) - < 

^A2/ 

0 

0 

mAs2 


^CA3 

J 


^A 
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[9) 


In order to simplify the terms in the following development, define 

,2 


[a] = 


IaS^ 
-Ha03S 
0 


Ha03S 

Ias2 

0 


0 

0 

mAS' 


( 10 ) 


and 


[t] = 


(/V2)rjn 

0 


-(V2)r„, 

-n/2)rn, 


(11 ) 
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Equation (9) becomes 




( 12 ) 


The bearing forces are functions of the bearing gap displacements {gx^ and 
commands {gxc^ as given by 


{^x} - Gf({9Xc} “ {^x}) 


( 13 ) 


where Gp is defined by equation (8) . Converting gaps to angles and displace- 
ments by using equation (1 ) results in 


Ti-y/jAcl _ r«A.' 




( 14 ) 


or 


([a] + Gp[T]'^[T] )|^~-^= Gp[T]'^[T]|j 


^AC_ 

*^CA3C 


(1 5) 


The system characteristic equation becomes 


det|[A] + Gp[T]'^[T]| = 0 


( 16 ) 


Expanding the determinant and making appropriate substitutions results in 

QmAs2 + 3 (KrS + Ka)]{[ias2 + ( 3 / 2 ) r„,2(KRS + Ka)]^ + Haq3S^} =0 (17) 

Making the thin rim assumption 


2 

= — 


( 18 ) 


I 



and setting ^^03 to zero allows equation (1 7) to be simplified to 


's2 + 


I"a/3 


s + 


.*)■ ■ • 


( 19 ) 


which illustrates the zero-mcanentum decoupling mentioned earlier. Referring 
to equation (1 7 ) , it can be seen that the translational roots are obtained from 

m^^s^ + 3 (KrS + = 0 


and are unaffected by the system momentum. The natural frequency of a single 
station becomes 



with a damping ratio of 


Kr 

2 ^KA(niA/3) 


In order to examine the effects of momentum, the second product term of equa- 
tion (1 7) is expanded to yield 


,4 ^ 3r^ ^3 ^ 

Ia 


Ha 03 3rm^KA 
Ia2 Ia 



s 


2 



KrKaS 


3 rn,2KA\^ 

2 Ia J 


0 


( 20 ) 


With the use of equation (18), equation (20) can be expressed in terms of the 
zero-momentum natural frequency and damping ratio of a' single station as 
follows: 
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4p(0nS- 


(^) 2p2)(On2 


+ 4P(j0j^3s + = 0 


( 21 ) 


Using the parameters of the laboratory model AMCD (described in the appendix) 
with set at the laboratory model design goal (ref. 5) of 262.7 N/iran 

(1500 Ib/in.) and assuming rated momentum (4068 N-m-s or 3000 Ib-ft-sec) , the 
damping ratio of the dominant system roots can be plotted against the damping 
ratio of a single bearing station at zero momentum as shown in figure 2. The 
implication of figure 2 is that if no practical limitations on rate gain exist 
then acceptable system performance is possible by designing for independent 
station control. 

The roots of equation (21), which represent the system poles, have some 
interesting characteristics which warrant discussion. At zero momentum there 
are, of course, two identical sets of poles- As the rim momentum increases, 
the poles form a set of high-frequency and a set of low-frequency complex poles. 
One set of poles increases in natural frequency as the rim momentum increases, 
and the other set decreases in natural frequency- The damping ratio is the 
same for each set of complex poles. For purposes of illustration, consider a 
damping ratio p of 1.65 which, from figure 2, should result in Pg of 0.7 
at full speed- With this value of p and the same parameters used to generate 
figure 2, the zero-speed poles of equation (21 ) become 


(s + 63.2)2(s + 554.9)2 


Figure 3 shows the motion of the low-frequency poles in the upper left quadrant 
of the s-plane as rim speed is varied from 0 to 2700 rpm- Figure 4 is a similar 
plot for the high-frequency poles. If the poles are complex at zero rim speed, 
both sets start at the same location in the s-plane and separate as rim speed 
increases. 


Radial Design 

The radial design problem is somewhat simpler than the axial design since 
momentum coupling is not involved. One consideration for the radial system, 
however, is rim growth at high speeds. A method of compensating for this in 
the laboratory model AMCD is presented in reference 5. As was the case with 
the axial system, the radial system will only be discussed in general terms. 

The equations of motion for the rim in the radial direction are 


\pr2J 
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where Frt and Fr 2 defined by equation (4) . In order to simplify the 

terms in the following development/ define 


[m] = 


1/2 

f^/2 

1/2 

-f/2 

-1 

0 

1 for 

a given : 


( 23 ) 


The force command for a given radial bearing F^(] can be written as 
^RC ^Fr(9rc “ 9r) 


where Gp^ is the forward loop transfer function of the radial loop/ gpQ is 
the radial gap command/ and gp is the radial gap. The force produced by a 
given bearing becomes 


^R = ^R^RC 


(25) 


where Ggp is the radial bearing transfer function. The total bearing forces 
become 


r 

^Ra 

^Rb 


A ^ 

^RaC ■ 9Ra 


) “ GprGbp /gRbC “ 9Rb > 


'Rc 


9RcC " 9Rc 


V. 
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The radial forces along axes 1 and 2 become (from eq. (4)) 




r 

9RaC ■ 9Ra 


) - ^FR^R^W^'^ \9RbC ■ 9Rb / 

I'saJ 


9RcC 


9Rc 


J 


(26) 


(27) 


1 0 



In terms of rim displacements ^ equation (27) becomes (from eq. (2)) 


r^Rll T , j’^CAlC ~ *^CA1 

) * GprGbr[m] [m]' 

|^CA2C ~ *^CA2J 




( 28 ) 


Using equation (22) , equation (28) can be written as 


(niAS^[l] + GprGbr[m] [m])| \ = GfrGbrLm] [m] 


rCA2, 


r ^CAl cl 

\5CA2C j 


(29) 


where [l] is the identity matrix. The system characteristics equation becomes 


detj(mAs2[l] + GfrGbr[m]^[m]) = 0 


(30) 


Expansion of this determinant results in 


^m^^s^ + 1 - 0 


(31) 


If perfect actuators are assumed (i.e., Gqr = 1) and if is a function of 

position and rate, then the radial system design reduces to that of a second- 
order spring mass system. 


PRACTICAL CONSIDERATIONS 

Initial tests with the laboratory model AMCD have provided insight into 
limitations of two hardware design approaches taken. These approaches, which 
are also discussed in reference 6, are (1 ) completely unidirectional layup of 
composite materials in the rim and (2) permanent magnet flux biasing for the 
bearings. Other hardware characteristics that warrant discussion are dis- 
tributed bearing elements and magnetic bearing power driver and power supply 
current limits. 


Rim Fabrication 

The laboratory model AM::d rim was fabricated by wrapping a graphite filament 
tape, impregnated with resin, on a special spindle producing 1 00-percent circum- 
ferential fibers (ref. 5) . As discussed in reference 6, this makes the rim 
susceptible to creep in a direction parallel to the spin axis; this results in 
deviations of the rim out of the spin plane which produces an equivalent dis- 
turbance input to the axial bearings as the rim spins. In the laboratory model 
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AMCD, Static deviations of the rim have been measured to be a maximum of 
1.651 mm (0.065 in.) to a minimum of 0.2032 mm (0.008 in.) depending on how the 
rim is stored. These numbers are total deviations (measured from peak to peak) 
and were measured midway between suspension stations with the rim suspended. 

The rim deviations, or warp, include two high and two low points per rim revolu- 
tion which approximate a sine wave disturbance at twice the wheel spin frequency. 
The origin of the basic two cycle per revolution shape of the warp is unknown. 

By experimenting with different methods of storing the rim, it was found that 
rim deviations could be held consistently to within 0.381 mm (0.015 in.) peak 
to peak. 


Permanent Magnet Flux Biasing 

As mentioned earlier, the AMCD laboratory model magnetic bearings utilize 
permanent magnet flux biasing. Advantages of this technique include: (1) A 

linear relationship between force and current at a given operating point can 
be easily obtained; (2) with permanent magnets supplying a ^rtion of the flux, 
the power required to suspend the rim in a 1 g (Ig = 9.8 m/s^ (32.2 ft/sec ^) ) 
environment is reduced. This technique does, however, present some problems 
from a control system standpoint. In order to discuss these problems, the 
linearized mathematical model of the AMCD laboratory model magnetic bearing 
actuator, derived in reference 9 (eq. (43)), is incorporated into equation (13) 
to yield 


{Fx} - KBGp({gxc^ ” 


(32) 


where the force command is represented as the current to the electromagnet. 
Equation (14) then becomes 






(33) 


and the system characteristic equation becomes 


det|[A] + (KbGf - %) j = 0 


(34) 


Expanding the determinant and making appropriate substitutions results in 


[mAs2 + 3 KbKrS + 3(KbKa - Km)]{[lAs2 + 3/2 r„2KgKRS 


(3/2)r„,2(K3KA - + ^ 03 ^^}= 0 


+ 


(35) 



Expanding the rotational part results in 


4 . 3rm2KBKRs3 

3 + — T" 


Ha 03 ^ 3rn,2 (kbKa - Q ^ U 




Ia 


2 Ia 


2 \2 

i *'bKr(KbKa - Km)s + 


3 _^ (kbKa - Km) 

Ia 


2 ""ni 


= 0 (36) 


Making the substitutions 


- _ J C^bKa - %) 

mp^/3 


0) 


n 


and 


KrK 


P = 


B^R 


2/{KbKa - Kro) (niA/3) 


equation (36) can be put in the following form: 


+ 4pfja^s^ + 


+ 2(1 + 2p2)SiJn2 


s2 + 4p/jaj^^s + = 0 


(37) 


which is the same form as equation (21 ) . As illustrated by figure 2, if no 
practical limit on rate gain exists, then acceptable system performance at full 
rim speed can be obtained. In the case of the AMCD laboratory model, an upper 
limit was imposed by a rim bending mode being driven unstable (ref. 5) . The 
upper limit on rate gain, from reference 5, was calculated to be 66.3 A-’S/m 
(20.2 A-sec/ft) with a position gain of 19571 A/m (5965 A/ft). This results 
in p = 0.365 (using the value of K^i from ref. (9)) and o)j^ = 1 57.86 rad/s. 
Using the relationship H;^q 3 = ^a34^s' curve of figure 5, which is a plot 

of the damping ratio of the dominant system poles versus rim speed, can be 
generated. As mentioned earlier in the discussion of equation (21), equa- 
tion (37) has two sets of complex poles which are equal at zero rim speed. As 
the rim spins up, the natural frequency of one set of poles becomes higher than 
the zero-speed value and for the other set it becomes lower. Figure 6 is a plot 
of the damped natural frequency of the poles of equation (37) versus the rim 
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spin speed. The damping ratio of both sets of poles is the same. The damped 
natural frequency, defined as 




( 38 ) 


is the frequency contributed to total system response by each set of poles and 
is the frequency at which the gain peaks for a set of underdamped poles on a 
magnitude frequency response plot. (For a discussion of frequency response 
methods and transient analysis by transfer function, see ref. 10.) Note from 
the figure that the low-frequency, or dominant, poles have a damped natural 
frequency a)j)L equal to twice the rim speed at approximately 500 rpra and equal 
to the rim speed at approximately 830 rpm. The crossover at twice rim speed 
is significant because the fundamental frequency of the rim warp is at this 
speed. Because of the increased motion of the rim in the bearing gaps at this 
resonance point, the linearization of the magnetic bearing actuators is no 
longer valid. Some insight, from a linear analysis standpoint, into the 
dynamics of large motions of the rim in the bearing gaps can be gained by refer- 
ring to figures 7 and 8 . These are plots of K 3 and versus their corre- 

sponding operating points in the bearing gaps. With rate and position gains 
fixed at 66.3 A-s/m (20.2 A-sec/ft) and 19571 A/m (5965 A/ft), and p vary 

over the same range of operating points as shown in figures 9 and 10. 


Distributed Bearing Elements 

As described in the appendix, each magnetic bearing station consists of 
four magnetic elements top and bottom for the axial suspension and four elements 
mounted around the inside of the rim for the radial suspension. The axial posi- 
tion sensor is mounted in the center of the top four bearing elements, and the 
radial position sensor is mounted in the center of the radial elements. This 
means that the bearing forces are being applied at points displaced from the 
point at which position is sensed. Incorporating the bearing element and sensor 
positions into equation (32) results in a negligible change in the system charac- 
teristic equations. The effects of rim warp and nonlinear bearing characteris- 
tics combined with distributed bearing elements produced a dynamic system model 
which was best investigated using the digital computer simulation described in 
reference 7. 


Power Driver and Power Supply Current Limits 

The laboratory model AMCD magnetic bearing power drivers (described in 
ref. 5) have a voltage limit of ±20 V. This voltage limit, in conjunction with 
the bearing coil resistance, results in an equivalent current limit of 13.13 A. 
In addition, the system power supply has a current limit of 30 A. 
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SIMULATION RESULTS AND DISCUSSION 


In order to evaluate the effects of rim warp, permanent magnet flux biasing, 
distributed bearing elements, and power driver and power supply current limits, 
the digital computer simulation, described in reference 8, was developed. Fig- 
ure n contains plots of peak-to-peak rim motion in a given axial magnetic bear- 
ing (peak-to-peak motion was the same in all stations) versus rim spin speed 
for the full nonlinear model with two rim warp amplitudes. The two amplitudes 
were 0.9525 mm (0.00375 in.) and 0.1905 mm (0.0075 in.) which produced peak-to- 
peak static deviations measured midway between suspension stations of 0.381 mm 
(0.015 in.) and 0.762 mm (0.030 in.), respectively. Figure 11 also contains 
plots of peak system current for the same rim warp amplitudes versus rim spin 
speed. For comparison purposes, simulation runs were made using a single point 
nonlinear axial bearing model and a single point linear axial bearing model, 
each with power driver and power supply current limits removed. The results 
are presented in figures 12 and 13. 

The main points to be made about these results are: (1 ) Rim motion is 

stable through the resonance excited by the rim warp for all cases. (2) The 
rim speed at which resonance occurs, for all cases, is very close to that pre- 
dicted by figure 6. (3) The resonance peak, for the rate gain considered, was 

approximately twice the warp peak amplitude. (4) The amplitudes and shapes of 
the plots for the single point linear and single point nonlinear bearings were 
very similar, and these plots showed only slight differences from the plots for 
the nonlinear distributed case. 

All of this tends to lend considerable confidence in subsequent linear 
analyses using the single point linear bearing model. However, the simulation 
runs differed from the actual hardware performance in two ways. First, the 
laboratory model was unable to attain a stable spin speed in excess of 475 rpm. 
Second, the radial current increase, as the rim was spun up, was much lower in 
the simulation than in the laboratory model. Recent tests, with a new magnetic 
bearing system installed in the laboratory model, have identified a static 
imbalance of the rim in the radial direction. Since the radial rim warp was 
the only anomaly included in the simulation, it appears likely that the major 
increase in current, as rim speed was increased in the laboratory model tests, 
was due to the static imbalance. 

In an effort to identify possible mechanisms for dynamic instability, simu- 
lation runs were made using the full nonlinear model with increasing values of 
warp amplitude. The system became unstable at a warp amplitude of 0.3226 mm 
(0.0127 in.). The speed at which the instability occurred was approximately 
770 rpm. The cause of the instability was a reduction in effective gains due 
to power supply and power driver limits being exceeded. It is interesting to 
note that the system became dynamically unstable before rim motion became large 
enough to trigger the backup bearing system. 

In order to investigate the effects of zero-speed single station damping 
ratios on system stability, further runs were made using the warp amplitude of 
0.3226 mm (0.0127 in.) but with zero-speed single station damping ratios of 
0.2 and 0.7. With a damping ratio of 0.2, the speed at which the system went 
unstable was reduced to approximately 410 rpan. When the damping ratio was 



Increased to 0.7^ the system was stable over the same range even though the 
power driver and power supply current limits were exceeded. 

To summarize, the results obtained from the simulation were: (1) With rim 

warp as the only anomaly included in the radial system, no power supply or power 
driver limits were exceeded for axial rim warp values up to 0,1905 mm (0,0075 in.) 
peak to peak, and the system exhibited stable operation. Results for distributed 
nonlinear, single point nonlinear, and single point linear axial actuators were 
similar. (2) By increasing the value of the axial warp amplitude, a mechanism 
for dynamic instability for the axial system was identified. This mechanism 
was the reduction in effective gains due to power supply and power driver limits 
being exceeded, (3) Dynamic instability occurred before rim motion became large 
enough to trigger the backup bearing system. (4) Under conditions where dynamic 
instability was possible, zero-speed single station damping ratio had a sig- 
nificant effect on the speed at which the instability occurred. 


CONCLUDING REMARKS 

A linear analysis and the results of a nonlinear simulation of a magnetic 
bearing suspension system which uses permanent magnet flux biasing has been 
presented. The magnetic bearing suspension is part of a 4068 N-m-s 
(3000 Ib-ft-sec) laboratory model AMCD, Included in the simulation were rigid 
body rim dynamics, linear and nonlinear axial actuators, linear radial actuators, 
axial and radial rim warp, and power supply and power driver current limits. 

The linear analysis of the magnetic bearing suspension system for the 
laboratory model AMCD indicated that stable system operation can be achieved 
by utilizing independent single station axial control. The single station 
bearing control system design can be performed using zero rim speed parameters. 
System performance at rated rim speed is influenced by the zero-speed single 
station damping ratio. 

The nonlinear simulation indicated that for rim axial warp amplitudes up 
to 0.1905 mm (0.0075 in.) peak to peak and for system operation where power supply 
and power driver limits were not exceeded, the linear and nonlinear axial bearing 
models gave similar results and the system was stable, as predicted by the linear 
analysis. For this warp amplitude range, however, the simulation runs differed 
from the actual hardware performance in two ways. First, the laboratory model 
was unable to attain a stable spin speed in excess of 475 rpm. Second, the 
radial current increase, as the rim was spun up, was much lower in the simu- 
lation than in the laboratory model. Recent tests, with a new magnetic bearing 
system installed in the laboratory model, have identified a static imbalance 
of the rim in the radial direction. Since the radial rim warp was the only 
anomaly included in the simulation, it appears likely that the major increase 
in radial current, in the laboratory tests, was due to static imbalance. 

By increasing the value of the axial rim warp in the simulation, a mech- 
anism for dynamic instability was identified. This mechanism was the reduction 
in effective axial system gains due to power supply and power driver limits 
being exceeded. Dynamic instability occurred before rim motion reached a level 
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that would have triggered the backup bearing system. Under conditions where 
dynamic instability was possible, zero-speed single station damping ratio was 
found to have a significant effect on the speed at which instability occurred 
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APPENDIX 


LABORATORY MODEL AMCD 

This appendix presents a brief description of the laboratory model AMCD 
assembly. A more detailed description of the subsystems are given in refer- 
ence 5. The laboratory model, shown in figure AT , consists of a graphite-epoxy 
composite rim which is 1 .6 m (63 in.) in diameter, weighs 22.5 kg (49.5 lb), and 
is designed to rotate at a speed of 2741 rpm. At this speed the rim momentum 
is 4068 N-m-s (3000 Ib-ft-sec) . The rim is suspended by three equally spaced 
support stations. Magnetic bearing elements located in the support stations 
interact with a low-loss ferrite material, embedded in the rim, to produce 
radial and axial suspension forces. Electromagnetic stator elements, also 
located at the support stations, push and pull against 72 equally spaced samar- 
ium cobalt permanent magnets, embedded in outer edge of the rim, to produce 
spin torques. The stator element drive electronics are commutated by signals 
from a Hall effect device which senses the position of the magnets. 

In order to prevent damage to the rim in the event of a magnetic suspen- 
sion failure during spin tests, the AMCD laboratory model includes a backup 
bearing system. The backup system includes six bearings (two per suspension 
station) which are designed to slow and support the rim. The bearing design 
provides hydrostatic air pads for radial control and hydrodynamic air pads for 
axial control. The backup and suspension bearing assemblies are attached to 
an aluminum baseplate. A vacuum cover (not shown) fits over the bearing-motor- 
rim assembly and also attaches to the baseplate. The cover is used for high- 
speed spin tests only. 

Figure A2 presents a close-up of one of the suspension stations. There 
are four magnetic bearing elements, top and bottom, for the axial suspension 
and four elements mounted around the inside of the rim for the radial suspen- 
sion. The center structure mounts the axial and radial position sensors for 
this station. The sensors consist of a variable impedance bridge with an 
active and a reference coil. They are sensitive to anything that changes the 
inductance of the active coil in the sensor probe such as the close proximity 
of a magnetic material. In the AMCD rim the ferrite material is sensed. Fig- 
ure A3 is a cross-section drawing taken through a bearing station and shows the 
rim, bearings, and spin motor elements in more detail. The magnetic bearing 
gaps with the rim centered are 2.54 mm (0.1 in.). Pertinent physical parameters 
of the laboratory model AMCD are listed in table A1 . 
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TABLE A1 

PHYSICAL PARAMETERS OF LABORATORY MODEL AMCD 


Spin inertia, I^^ 3 , kg-m^ (slug-ft^) 14.36 (10.6) 

Transverse inertia, I/^, kg-m^ (slug-ft^) 7.18 (5.3) 

Mean radius of rim, rj^, m (ft) 0*8 (2.625) 

Mass of rim, m;^, kg (slugs) 22.48 (1 .54) 


Figure Al 


1 . 
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laboratory model. 
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Figure 6*- Variation of damped natural frequency of system poles 

with rim speed. 
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